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We study a current transport through a system of a few grains connected with tunneling links. 
The exact solution is given for an arbitrarily connected double-grain system with a shared gate in 
the framework of the orthodox model. The obtained result is generalized for multigrain systems 
with strongly different tunneling resistances. We analyse the large-scale nonlinear conductance and 
demonstrate how the sample topology can be unambiguously deduced from the spectroscopy pattern 
(differential conductance versus gate-bias plot). We present experimental data for a multigrain 
sample and reconstruct the sample topology. A simple selection rule is formulated to distinguish 
samples with spectral patterns free from spurious disturbance caused by recharging of some grains 
nearby. As an example, we demonstrate experimental data with additional peaks in the spectroscopy 
pattern, which can not be attributed to coupling to additional grains. The described approach can 
be used to judge the sample topology when it is not guaranteed by fabrication and direct imaging 
is not possible. 

PACS numbers: 73.23.Hk, 73.50.-h, 85.30.Vw 
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I. INTRODUCTION 



The electrical properties of nanostructures with ever-smaller size have been studied intensively during the past 
decade in a search for new quantum effects and novel types of electronic devices. The scientific interest has moved 
towards objects with sizes of a few nanometers. Such objects cannot be prepared by conventional electron beam 
lithography with a resolution being limited to ~ 20 nm. The laboratory methods, developed thus far, make use of pre- 
fabricated nanoscale objects: semiconductor jO or metallic nanocrystals made by colloidal chemistry, organometallic 
carborane clusters J3, and the fullerene molecules The granular films formed due to surface atom diffusion at 
room temperature or in a process of an avalanche reconstruction of a quench-condensed precursor layer JTot were 
also used as a natural source of nanoparticles. These prefabricated nanoparticles have been spread on a surface and 
studied with STM [|| ; or placed into a gap between two macroscopic electrodes defined by electron beam litography Q| 
or made by a break-junction technique Q. 

The final sample geometry, i.e. the number of nanoparticles in a gap and their arrangement is not under experimental 
control. With some luck, one can get a sample where a single nanograin bridges two macroscopic electrodes. Typical 
yields of 25-30% have been reported ||]8) . 

Several methods to get a better defined sample geometry were suggested. A hole in a silicon-nitride membrane 
was used to restrict a contact area [ pH7| . A self- assembling of gold clusters on electrodes covered with monolayers of 
organic molecules was used by Ref. Q2J. An alternative way is to modify a sample by in situ material deposition, as 
described in |s|Jlot|. 

Nevertheless, not a single method ensures, just by the fabrication procedure itself, that there is no more than one 
nanoparticle in a gap. Moreover, even if the current mainly goes through a single grain, there may be some additional 
grains nearby which can potentially affect transport properties due to capacitance coupling. Sometimes, the sample 
geometry can be imaged directly by STM or SEM, like in Jl],||], but sometimes not, as the samples are metastable and 
can only exist at cryogenic temperatures fl(i|j , or nanoobjects are simply too small, like the fullerene molecules [Q. 

The only universal way to judge the sample topology is to analyze current transport data. Several criteria have 
been suggested in the literature, all based on the Coulomb blockade phenomena |fj"lfl . A nano-scale object, coupled to 
macroscopic source and drain electrodes, is a Single Electron Transistor - SET (although not every method allows to 
implement a third electrode - gate). The exact periodicity of the conductance variation versus bias and gate voltages 
in a single grain SET is violated if more than one grain is involved in the current transport [ fl2"| . So one can formulate 
an intuitive "selection rule" for single grain samples as "the ones showing Coulomb-staircase structure" (^,§,0] or, more 
precisely, " the ones with a single set of equally spaced Coulomb-staircase peaks in the differential conductance" Q . 

Within the scope of the orthodox model for SET the Coulomb staircases are the only peculiarities in the 
differential conductance. The real nanoscale samples, however, never follow this model exactly. First, if a typical 
level spacing S between the quantised energy levels in a nanoparticle is comparable with the charging energy Eq, 
then the Coulomb diamonds have different sizes and the Coulomb staircases are no more periodic |0|. Second, if the 
temperature is below the inter-level spacing 5, then every discrete energy level in a grain gives an additional peak in 
the differential conductance dl/dV Q. In fact, it was the discrete excitation spectrum in a nanoparticle, which was 
the primary subject of the research in references |^ S|. 

Therefore, one should distinguish between the peculiarities coming from the discrete energy levels and a potential 
disturbance due to additional grains nearby |l4| . If all energy levels in the grain energy spectrum are equally shifted 
by the gate potential then any discrete energy state is revealed in a dI/dV(V, V g ) plot (V is the gate and V g is the 
bias voltage) as a straight line parallel to the diamond edge [||. This allows us to reformulate the selection rule for 
single grain samples as "the ones with all spectral lines in a dI/dV(V, V g ) plot parallel to the diamond edges". 

However, the gate voltage shifts all discrete levels the same way if only the external electrostatic field is effectively 
screened within a thin surface layer A <c R, where R is characteristic sample radius. This is valid in the extreme case of 
a macroscopic metallic sample, but definitely does not take place for the opposite extreme case - a fullerene molecule. 
Moreover, additional spectral features can be expected due to excitation of mechanical degrees of freedom PUDs!, 
including the shuttle effects predicted in |l6f| , and they could have a distinct gate dependency. To summarize, no 
single qualitative criterion, suggested thus far, provides a conclusive argument to prove a single grain sample topology. 

In this paper we report the exact solution (within an orthodox model) for an arbitrary double grain system, and 
some mathematical results concerning the multigrain systems. Based on our analysis, we can formulate an excluding 
criterion, i.e. we can specify what kind of spectral pattern can never be the result of spurious disturbance from one 
or more additional grains nearby. We also present, as an example, experimental data of this type. 
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II. ENERGY DECREMENTS ON TUNNELING 



In our analysis we follow the solution of an orthodox model for a single-grain SET, as given in |TT[ , with a proper 
generalisation for a double-grain system. Whenever possible, we also use the notations introduced in Q. 

The electrostatic charging energy decrement AK a jj during the tunneling jump of an electron between the grain a 
and the lead electrode (3 can be found directly and equals (see also Fig. [j] for notations): 
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Here we defined Cei = C n + C 12 + C g \ + C, Ce2 = C 2 i + C22 + C 92 + C, and = CeiCe2 - C 2 . We also defined 
a = 2 if a = 1 and a = 1 if a = 2; the same definition for (3. 



Equations (2.1) and (2.2) are highly symmetric - see Appendix A. 



III. QUASIPERIODICITY 



The position of Coulomb staircases in a (V^ V^) plot is determined by the condition that the electrostatic energy 
decrement on tunneling through some barrier crosses zero: 
A^ a/3 (m,n 2 ,y,K) -0. 
From (2.1) and (2.2) we have: 
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Because of time-reversal symmetry between tunneling jumps a — > j3 and /3 — ► a, the conditions AS* Q ^(n Q ) = and 
A^ Q £j(n Q + 1) = impose the same relationship on V, V g , and only 5 distinct sets of pecularities are defined by (|3.l| ). 
We can draw some important conclusions from (3.1): 



a. All Coulomb pecularities form straight lines in a gy (V, V^) plot. Indeed, any equation in (3.1) has a linear form 



A a pV + BapVg — CafjUi + D a pn 2 + const. (3.2) 

b. Each tunneling barrier defines a set of pecularities which are parallel to each other, i. e. a double dot system 
has maximum five distinct slopes in the ^ pattern. 

c. There is no more periodicity with respect to a gate voltage, like in a single-dot system. Nevertheless, the whole 
pattern is quasiperiodic, because all staircases coming from the first grain have offsets 1%% + (C / C-£ 2 )n 2 + • • •, from 

the second one - (C/Csi)^i + 7i 2 + . . ., and from intergrain tunneling - (Cs2 — C)ni — (Csi — C)7i 2 + • • •• 

IV. MASTER EQUATION FOR EVOLUTION 

If both grains have a continuos spectrum, then for tunneling rates r Q ^(ni, n 2 ) through the lead-grain gaps we have: 
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and for intergrain tunneling rates: 
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where /^(V) is the current- voltage characteristics of the corresponding tunneling barrier. Assuming Ohms law, we 
put I a fj{V) = G a pV, but all formulas can be easily generalized for the case of a discrete energy spectrum in the 
grains. 

The master equation for the double dot system is 
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And the current is given by 
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V. NUMERIC SOLUTION OF A DOUBLE DOT SET 



Master equation (4.3) can be solved numerically in a straightforward manner: for any given (V,V g ) one should 
trace the evolution of p n ,m until an equilibrium distribution is reached (see Appendix B for the details of the numeric 
procedure), the current is then given by (4.4). This way the current can be computed for an arbitrary double grain 
system. 

To test the program, we solved a general double dot system and compared the program output with the result of a 
Monte-Carlo simulation obtained with a commercial single electron circuit simulation software - SIMON 2.0 Jj"7| . As 
one can see, all Coulomb peculiarities have identical positions in both Fig. ||a and Fig. ||b. 



VI. DOUBLE DOT SYSTEM WITH STRONGLY ASYMMETRIC TUNNELING BARRIERS 

Fig. ^illustrates the common rule: in a multigrain system any tunneling gap gives a set of peculiarities on -^(V, V g ) 



plot with one and just one unique slope. Equation set (3.1), however, only specifies possible staircase positions, not 



amplitudes. The amplitudes depend on probability distribution p ni ,n 2 as defined by (4.3). 

For example, lets consider the case of negative bias voltage V with electrons flowing in Source-to-Drain direction. 
Whenever the small variation of V and/or V g makes, say, (n, to) positive, electron from Fermi level in a Source 

lead can tunnel into a Fermi level in the first grain. One additional intermediate charge state {n + 1, m} becomes 



accessible, resulting in one extra term in current sum (4^4) due to to {n,m} — > {n + 1,to} transitions. The resulting 



current increment, however, is proportional to the probability of a system to be found in a {n, to} state multiplied by 
{n, m} — > {n + 1, to} tunneling rate: T n(n, m)p nTn . If p n<m is small, i.e. the {n, to} state is statistically unfeasible, 
then a newly opened channel does not contribute much the total current, and the corresponding staircase is suppressed. 

For clarification, we can refer to a single grain SET. In a symmetric case (G\ = G2) all charge states, accessible 
from the energy considerations, have comparable feasibility: p nmax ~ Pn max ~i ~ • •• ~ P-n max - On the contrary, in 
the asymmetric case G\ 3> G2 the grain is almost in equilibrium with the low resistance electrode. In a rare event of 
tunneling through the low-conductive gap G2 equilibrium is violated, but almost immediately restored by tunneling 
through a high-conductive barrier G\. Two successive tunneling events through G2 not interleaved by tunneling 
through Gi are even less feasible. As a result, the charge fluctuations are strongly suppressed and in the G2 ^ 
limit the probability distribution reduces to a single charge state n max : p n — § n ,n ma x- 

Whenever a bias or gate increment makes tunneling through the G\ gap into the n+ 1 state possible, an equilibrium 
grain charge switches to e(n+ 1), and the current is incremented stepwise. On the contrary, when the voltage variation 
opens a new tunneling channel through the G 2 gap, the equilibrium charge remains the same, and no current step 
takes place. While a symmetric SET has two distinct staircase families, the first one defined by Ki(n, V, V g ) — and 
the second by ^(n, V, V g ) = 0, in an asymmetric system one of them is suppressed. 

Exactly the same kinetic effect takes place for a double grain (multigrain) SET: if tunneling gaps have strongly 
different conductances G a p, then the vast majority of accessible charge states becomes unfeasible, and for any given 
pair (V, V g ) just one charge state {rii(V, V g ), ri2(V, V g )} dominates. If only V and/or V g variation changes the dominant 
charge state, a prominent staircase appears. For any grain a, the grain charge is defined by the most conducting 
tunneling gap /3 max , and the charge state changes when K a p max (ni, n,2, V, V g ) crosses zero. All other staircases, 
corresponding K a ^p max are supressed. 

To summarize, in a multigrain system each tunneling gap defines one staircase family with a distinct slope in the 
(V, V g ) plane. But in a strongly asymmetric system only one staircase family per every grain survives - the one 
corresponding to the most conductive link. Although in the double grain system there may be up to five distinct 
staircase families, in the asymmetric case only two will stay - the first one corresponding to the n\ increments and 
the second one to ri2 increments. Three-grain system will have at least three staircase families and so on. And vice 
versa - a system with n distinct staircase slopes has no more than n grains. 

This qualitative consideration is confirmed by the numeric result presented in Fig. ||[ 



VII. A SYSTEM WITH TWO STAIRCASE FAMILIES: POSSIBLE TOPOLOGIES 

As it was discussed above, a system with two staircase families may have maximum two grains. Each grain must 
have one strong tunneling link with conductance much higher than the other tunneling gaps adjacent to this grain. 
This allows us to list all possible topologies directly. We note first that each grain has three links, and any one could 



5 



be dominant. Therefore, for two grains we have 3x3 different possibilities, all shown in Fig [|. (note that for generality 
we assume that the grains have different size). Variants a, c, e, and g are obviously symmetric to b, d, f, and h 
respectively. Variant i has no more than two staircase families if only one of the grain-electrode gaps has much higher 
conductance than the others. In this case the system reduces either to d or f. To summarize, we have four distinct 
cases: b, d, f, and h; we shall consider each one in more detail belowQ. 

Grains are linked to different electrodes - the case in Fig. |]b. Two visible staircase families are those coming 
from the gaps Gn and G*22- The dl / dV (V, V g ) plot looks slightly different depending on which of the remaining 
tunneling gaps - G12, G21 or G has the highest conductivity. The current either flows mainly through the big grain 
(G21 > max(G, G12) - the numeric result is shown in Fig. ||bl), through the small one (G12 > max(G, G21) - see 
Fig. ||b2), or two grains connected in series (G > max(G2i,Gi2) - presented in Fig. ||b). Note that in all cases the 
Coulomb staircases have exactly the same position in the (V, V g ) plane, but the staircase amplitudes (and even the 
signs) are different. The regions of the Coulomb blockade are also different. In the general case of G12 ~ G21 ~ G 
the Coulomb blockade is lifted as soon as at least one of these current channels opens for transport. 

A big grain is linked to a small one, which is linked to an electrode - the case in Fig. ||d. The current either goes 
through the big grain (G22 > G12) or bypasses it (G22 < G12 ). The corresponding dl /dV plots are shown in Fig. ||dl 
and Fig. ||dl respectively. 

A small grain is linked to a big one, which is linked to an electrode - the case in Fig. ||f. Two possible topologies 
complementary to dl and d2 are presented in Fig. ||fl and Fig. [|f2 together with dl /dV plots. 

The last possibility - both grains are linked to the same electrode - is shown in Fig. ^h. If the current goes mainly 
through the small grain, we have the dl/dV pattern presented in Fig. ||hl, otherwise - the one in Fig |]h2. 

In any case dl/dV pattern has a simple structure: there are big areas shaded by parallel and equdistant spectral 
lines (one is shown by dotted line in Fig. |^hl). The reason for this highly regular structure is quite obvious for 
topologies in bl, b2, hi, h2 - a single grain SET coupled to a Single Electron Box (SEB). Indeed, until SEB charge 
Qseb remains constant, the SEB affects SET in essentially the same way as a constant background charge (of course, 
there is also some capacitance renormalisation) . So all SET staircases are periodic with (V, V g ). Whenever Qseb 
changes, the SET feels a different "background charge" and the phase of its Coulomb blockade oscillations is shifted. 
Remarkably, this high regularity survives even if two grains are connected in series, as in dl and fl. 



VIII. COMPARISON WITH EXPERIMENT 



The analysis presented in this paper was to a large extent motivated by our recent experiments with nanocon- 
strictions in Quench-Condenced (QC) films |9|,[l0|]. We have developed a method to define constrictions with width 
and length of about 5 nm, bridged to two macroscopic electrodes. At small thiknesses QC films are granular jl^], so 
there are typically a few metallic grains separated by tunneling gaps within such a constriction - see Fig. ^| for the 
sample geometry. By additional in situ film depositions we merged some clusters together, finally reducing the sample 
topology to a single grain SET. The only opportunity to follow the sample topology was the analysis of the sample's 
dI/dV(V, V g ) pattern. We used the technique presented here to determine the number of grains in the constriction, 
to reconstruct the sample topology, and to extract all grain capacitances. The fast computation time (compared to 
the Monte-Carlo method) makes this procedure practical. 

The strong exponential dependence of the tunneling conductance on the gap width greatly simplifies the spectral 
pattern. Indeed, we can assume that the tunneling gaps between the grains are randomly and uniformly distributed 
from 1 A to 10 A. Then the tunneling resistances are randomly and uniformly distributed in logarithmic scale from 
10 kQ and up to 1 GQ. This virtually eliminates the kinetic charge fluctuations, and washes out all but one staircase 
family per grain, i. e. every staircase family identifies the unique grain. 

A typical experimental data are shown in Fig. [7] together with a numeric fit. The reconstructed sample topology 
is presented in Fig. [?]c: two grains, with charging energies e 2 /2Gs about 9 and 34 meV, are connected in series; and 
there is a small single electron box with charging energy ~ 52 meV nearby. 

Note that the inter-grain capacitance G is always small compared to coupling to the leads or the gate electrode (see 
caption to Fig. uft. This is because an inter-grain interaction is effectively screened by the common gate electrode. For 



Note that it is natural to assume in b, d, and f cases that a smaller grain is closer to the left electrode, and the bigger grain 
- to the right one. That's why all numeric results presented below were computed for a system with Cn > C12 and C21 < 622- 
We put this constraint to make the staircase slopes more realistic, it does not affect the cardinal pattern of dl/dV plot. 
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the sample in Fig. fj] the gate oxide layer was about 3 nm thick, so for two grains bigger than ~ 3 nm we always have 
C < C g . Moreover, any two grains placed more than a few nanometers away virtually do not interfere. This allows to 
apply our analysis of a double-grain system recursively, and to identify all grains in a multi-grain system one by one. 

Whenever an additional deposition merges two grains, one staircase family disappears. This evolution eventually 
leads to a single grain SET, and after a few additional depositions, to a single tunneling gap separating two macroscopic 
electrodes JToj ] . 

Intriguingly, for samples with very high charging energies (> 100 meV), we have always observed some non-orthodox 
spectral lines like the ones marked with 7 in Fig. ^. All lines 7 have the same slope, different from the slope of the 
main Coulomb staircases. As it was shown previously a system with two distinct slopes on the dl /dV pattern has 
no more than two grains involved in current trunsport. Would the lines 7 be due to recharging of some other grains, 
they all must come from recharging of one and the same grain. In this case they must be equidistant, which obviously 
does not take place. 

IX. CONCLUSION 

We proved the following statements for a multigrain system in the framework of the ortodox model: 

1. All Coulomb peculiarities in differential conductance form straight lines on the gate- bias plane. 

2. Each tunneling gap defines a set of staircases, which are parallel to each other. In an asymmetric system with 
strongly different tunneling resistances only one distinct staircase family per every grain survives - the one 
corresponding to the most conductive link; all other staircases are suppressed. 

3. There may be a maximum of n grains in a system with n staircase slopes. 

We presented the exact solution for an arbitrarily connected double-grain system. We demonstrated that although 
the zero-bias conductance is stochastic JEJ, the large-scale spectroscopy pattern dI/dV(V, V g ) (differential conduc- 
tance as a function of bias and gate voltages) is regular - all Coulomb staircases are quasiperiodic. 

We listed all possible topologies for a system with two different staircase slopes and presented the corresponding 
spectral patterns. 

Finally, we proved the following criterion: 

If there are only two slopes a and /3 for all spectral lines on spectroscopy pattern, and the lines with the slope /3 are 
not equidistant within the regions between two adjacent lines alpha, then they can not be attributed to the presence 
of one or more additional grains in the sample^. 

To conclude, we have analysed a multigrain system with a shared gate and found that the sample topology can be 
unambiguously deduced from dI/dV(V, V g ) pattern. We demonstrated how to distinguish peculiarities in tunneling 
spectra of single-grain samples from possible disturbance caused by some extra grains. Presented analysis can be 
applied for numerous systems, where the experimentalist is not able to image the sample directly. 
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XI. APPENDIX A: SYMMETRY CONSIDERATIONS 



Equations (2.1) and (2.2) are invariant with respect to the following symmetry transformations. 



2 This criterion can be generalized in a straightforward manner to the case of n slopes on dl/dV pattern, but this situation 
seems to be of low practical importance. 
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a. We can swap indices marking the grain numbers, i. e. we can redeclare grain no 1 as 2 and vice versa. So as a 
result of simultaneous replacement 

n>i ^> ri2, Cxi ^ C21, C12 ^ C22, C g i *=> C S 2, and Cm C^2 
we will have the following transformations: 

AKiff(ni,n 2 ) *=> A~lt 2 p(n 2l A^i /3 (rii,n 2 ) ^ A < K 2 f3{n 2 ,ni), and 
AKi 2 (ni,n 2 ) ^ AK 2 i(n 2 ,ni). 

b. Left/right symmetry: we can flip the whole network horizontally with simultaneous inversion of the bias voltage, 
which is equivalent to the following substitution: 

C\\ C12, C21 ^ C22, ^ —V. 
As a result, the tunneling jump from, say, the left tunneling barrier will be transformed into a jump through the right 
barrier in the opposite direction, so as 

AJ? a i(V,vg *=> a5? q2 (-V,V 9 ), A^fV.Vj) ±5 A^ o2 (-V,V fl ). 

Energy decrements for intergrain tunneling Aif^ and Ai^i are not affected by this transformation and they will 
convert to themselves. 

c. Time inversion symmetry: if for any tunneling event from state n\,n 2 to state n'^n^ the energy decrement is 
AK, then the inverse tunneling process from n'i, n' 2 to m, n 2 must have the opposite energy decrement — AK, and so 
A^n(ni,n 2 ) = -A^n(rii + l,n 2 ), A# 12(711, n 2 ) = -AS\ 2 (ni - l,n 2 ), 

Alt 2 i(ni,n 2 ) = ~A*K 21 (ni,n 2 + 1), A^K 22(711, 7i 2 ) = -AK~ 22 (m, n 2 - 1), 
and AKi 2 (ni,n 2 ) = -AK 2 i(rii — l,n 2 + !)• 



Any of a-c symmetries reduces the number of independent formulas in (2.1) by a factor of 2, so there is really just 



one independent equation in (2.1) 



There is yet another symmetry in (2.1) and (2.2) 



d. All energy decrements change the sign and the tunneling direction on simultaneous inversion of all voltages and 
charges, i. e. 

AK a p(ni,n 2 ,V,V g ) = -AKa^-n^ -n 2 , - V, -V g ), and 
A < K al3 (ni,n 2 ,V,Vg) = -Alt a p(-ni,-n 2 ,-V,-Vg), 



If we apply this transformation to the expression for current (4.4) then the current will change sign 



J(V, V g ) = —I{—V, -V g ), and ^(V, V g ) = §jH-V, -V g ). ( I I.I ! 



So the I(V, V g ) plot is antisymmetric, and the g^(V, V g ) plot is symmetric. In the presence of a background charge, 
the relative phase of Coulomb oscillations on the first and the second grain will be shifted and this symmetry will be 
violated. 

If we complete transformation d with the inversion of the electron charge e — -> —e then the current will not change. 
This could be reformulated as follows: the current through the system does not depend on the sign of charge carriers, 
i. e. whether the current transport is carried by electrons or holes. 



XII. APPENDIX B: NUMERICAL RECIPIES 



The master equation (4.3) is a system of linear differential equations. If we present the probability distribution 
Pni,n 2 as a vector p = \pj) where j is an index running over all (m, n 2 ) pairs, then we can rewrite (|4.3|) as 



P 



-Tp. 



(12.1) 



The formal solution of (12.1) is 



m 



,-*r 



p(o) = (e - tr + — r 2 
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-r 5 



24 



r 4 + 



(E - tf )(E - If )(E - |f )(E - if )p(0) + o(t 4 ), 



(12.2) 



where E is the unity matrix. 



As one can see from (12.2), the standard fourth-order Runge-Kutta method reduces in our case to four recursive 

t t t 

4' 3' 2' 



calls to a first-order Eiler's extrapolation formula with steps 4, 4, 4, t 
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To minimise execution time, it is better to use an adaptive step t. Ideally, it should be inversely proportional to 
the maximum eigenvalue of the evolution matrix T X max . Practically, it is good enough to replace an unknown X ma x 
with the sum of all eigenvalues J^Aj — * r (r), where the trace of the evolution matrix can be found from (13) and 
equals 



tr(T) 



nl,n2=-oo a,/3=l,2 



(12.3) 
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FIG. 1. Top: system topology and energy decrements on tunneling jumps; Bottom: equivalent schematics. 
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FIG. 2. General case of a double dot SET: differential conductance plot as a function of bias and gate voltages, a.- The 
result of Monte-Carlo simulation, b.- Direct numeric integration of master equation. System parameters: Cn = 16 x 10 -20 F, 
C 12 = 8 x 1(T 20 F, C21 = 9 x 1CT 20 F, C22 = 13 X 1CT 20 F, C gX = 8 x 10" 20 F, C g2 = 11 x 1(T 20 F, C = 10 x 1(T 20 F, 
Ru = 3 X 10 9 R 12 = 2x 10 9 n, R 2 i = lx 10 9 ft, R22 = 2 x 10 9 n, R = 20x 10 9 n. 
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FIG. 3. a.- Double dot system with identical tunneling conductances: Gn = G12 ~ G21 ~ G22 = G. There are five distinct 
staircase slopes, b.- Every grain has one well-conducting tunneling gap which governs the grain charge: Gu = G22 = 30G; 
G12 = G21 = 0. All but two staircase families are suppressed. Same capacitors for symmetric and asymmetric systems: 
Cn = 45 x 1(T 20 F, C12 = 6 x 1(T 20 F, C21 = 18 x lO" 20 F, C22 = 126 x 10" 20 F, C gl = 19 x 10" 20 F, C g2 = 47 x 10" 20 F, 
C = 37 x 1(T 20 F. 



a b 




FIG. 4. All possible topologies for a double grain system with no more than two staircase families. The case a produces two 
families, associated with tunneling gaps G21 and G12, the case b - families associated with Gu and G22, and so on. 
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-IV IV 

FIG. 5. dI/dV(V,V g ) patterns for a double-dot SET: all possible cases (together with Fig. Hb) when just two staircase 
families are present. For any topology sketch one solid line indicates the most conducting link for the bigger grain and another 
solid line - the most conducting link for the smaller grain. The dotted line shows the main current channel. Note that in cases 
bl and b2 the second grain gives the staircases with a negative differential conductance. There are also small areas of negative 
dl/dV in d2, f2, hi, and h2. 
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FIG. 6. Sample geometry: Thick low resistive electrodes are bridged by 5 x 5 nm 2 network of metallic clusters. Beneath is 
an electrostatic gate electrode made of the oxidized Aluminium. 
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FIG. 7. A typical example of experimental data with the corresponding numeric fit. a.- Two grains connected in series. 
The widest tunneling gap is between the small grain and the electrode - the same topology as in Fig. gfl. Inset: a numeric 
fit with parameters G21 = G = 30Gi 2 ; Gu = G22 = 0; G11 = 8 X 1(T 20 F, C12 = 134 x 1(T 20 F, G21 = 380 x 10~ 20 F, 
C22 = 114 x 10" 20 F, Ggi = 53 x 1(T 20 F, C g2 = 352 x 1(T 20 F, C = 41 x 10~ 20 F. b.- The same sample as in a. Current 
gradient across the direction indicated by an arrow was computed to highlight the third staircase family almost unresolvable 
in a. The third grain is coupled the same way as in Fig. |Hb2. Inset: a numeric fit with parameters G11 = G22 = 3OG12; 



G21 = G = 0; 
C g2 = 65 x 10" 



Gn 



= 45 x 10~ z ' 
G = 7 x 10" 



F, G12 
20 R 



120 x 10" 
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F, G22 = 44 x 10" /u F, G, 
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Gate Voltage (V) 

FIG. 8. SET with a charging energy ~ 90 meV. Non-orthodox lines 7 cannot be attributed to the recharging of some other 
grains. 
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